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Abstract 

The Ricci flow is a heat equation for metrics, which has recently been used to study the topology 
of closed three manifolds. In this paper we apply Ricci flow techniques to general relativity. We 
view a three dimensional asymptotically flat Riemannian metric as a time symmetric initial data 
set for Einstein's equations. We study the evolution of the area A and Hawking mass TW// of a 
two dimensional closed surface under the Ricci flow. The physical relevance of our study derives 
from the fact that, in general relativity the area of apparent horizons is related to black hole 
entropy and the Hawking mass of an asymptotic round 2-sphere is the ADM energy. We begin 
by considering the special case of spherical symmetry to develop a physical feel for the geometric 
quantities involved. We then consider a general asymptotically flat Riemannian metric and derive 
an inequality -^A'^^'^ < — 247r^/^A4iy which relates the evolution of the area of a closed surface S 
to its Hawking mass. We suggest that there may be a maximum principle which governs the long 
term existence of the asymptotically flat Ricci flow. 

PACS numbers: 02.40.-k,04.70.Dy 
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I. INTRODUCTION 

The Ricci flow 

mm 

has been used by mathematicians to understand the topology 
of three manifolds. It appears likely that these mathematical developments will also be 
useful in physics in the study of geometric theories like general relativity. Several papers 
have appeared dealing with physical applications of the Ricci flow. The 
Ricci flow (RF) is a (degenerate) parabolic differential equation, and is very similar to the 
heat equation. In a previous paper [l^ we explored an analogy between the Ricci flow and 
thermodynamics. The analogy is based on the observation that the Ricci flow (like the heat 
equation) loses memory of initial conditions, just as a physical system loses memory of its 
initial state as it approaches thermal equilibrium. As was noticed there, a slight modification 
of the Ricci flow yields Schwarzschild space as a fixed point. In this paper we look at the 
unmodified Ricci flow to see how some physically interesting quantities evolve with the flow. 
Energy and entropy are quantities of physical interest from the thermodynamic point of 
view. In general relativity these quantities take on a purely differential geometric meaning: 
the entropy is related to the area of black hole horizons and the energy to the ADM mass at 
infinity. We investigate the the evolution of these quantities under the Ricci flow and derive 
some inequalities relating them. 



The subject of entropy bounds 
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14l | has been of wide and sometimes controversial 



|15l | interest to physicists. It is believed that such bounds contain clues towards quantum 
gravity. Entropy bounds motivate ideas like "the holographic hypothesis". There is a 
geometric entropy bound proposed by Penrose, which connects two geometrical quantities: 
the area of the outermost horizon and the ADM mass. In the seventies, Roger Penrose 16| , 
in an attempt to test the idea of cosmic censorship (which is part of the "establishment 
view" of gravitational collapse), used physical arguments to deduce from cosmic censorship 
an inequality relating the ADM mass of an initial data set for GR and the area of its 



outermost apparent horizon: Madm > vIGttA. This inequality, which is saturated by the 
Schwarzschild space, has a clear thermodynamic interpretation. It states that Schwarzschild 
spacetime maximises geometric entropy for a given energy. Maximising entropy for fixed 
energy is a property which characterises a thermal state in statistical mechanics. The thermal 
character of Hawking radiation from a Schwarzschild black hole is entirely consistent with 
this interpretation. Penrose's inequality appears to capture something deep about general 
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relativity, with a thermodynamic statement expressed in geometric terms. In this respect, it 
is similar to the the area theorem ~v\ . which represents the second law of thermodynamics 
in geometric form. 

A counterexample to Penrose's inequality would imply a flaw in the establishment view. 
No counterexample has so far been found. Special cases of Penrose's inequality have been 



proved 1 18l. Il9l . l20l | using geometric flow techniques. Jang and Wald [21] showed that if one 



assumes the existence of the inverse mean curvature (IMC) flow, (a particular diffeomor- 
phism of a spatial slice), one could address Penrose's inequality in the special case of time 
symmetric initial data. Their work was based on the monotonicity of the Haw king mass 
under the inverse mean curvature flow, which follows from earlier work by Geroch 2^. This 
suggests a proof of the Penrose inequality [21]. The only gap in this proof is the existence 
of the IMC flow. This gap was fllled by Huisken and Ilmanen [2(3] who showed that the 
Hawking mass remained monotonic under a discontinuous version of the IMC flow. 

The motivation behind the present work is to explore whether the Ricci flow could lead to 
a new proof of the time symmetric Penrose inequality. The general Penrose inequality is still 
an open conjecture and new lines of attack (even in the time symmetric case) are certainly 
of interest. Is it possible that a smooth combination of the Ricci flow plus diffeomorphisms 
exists so that the evolution of Hawking mass is monotonic? The hope underlying the use 
of the RF is that the discontinuities (like those of the IMC flow) will be smoothed out by 
the RF term. There are known examples in fluid mechanics where discontinuities of inviscid 
flow are smoothed out by viscous effects. With this motivation of developing a Ricci flow 
approach to the Penrose inequality, we have made a beginning in this paper by studying the 
evolution of area and Hawking mass under the Ricci flow. 

In section II we give a brief review of the Ricci flow. For a more complete and rigorous 
account the reader is referred to the mathematical literature 23|, |2J] . Section HI lists the 



geometrical quantities of physical interest in this paper. Section IV paves the way for the 
general treament by treating the special case of spherical symmetry. This special case is 
useful since one can explicitly work out geometrical quantities of interest and develop a 
physical feel for them. Spherical symmetry is a good source of physical examples and coun- 
terexamples which guide the general study. In Section V we give up spherical symmetry and 
treat the evolution of area of a closed surface under the Ricci flow and derive an inequality 
relating the rate of change of area and the Hawking mass. Section VI treats the evolution of 
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Hawking mass under Ricci flow and puts forward a conjectured maximum principle which 
governs the long time behaviour of the flow. Section VII is a concluding discussion. Our 
metric conventions are from Ref. 25 1. 



II. THE RICCI FLOW 

Let (S, hab) be an asymptotically flat, three dimensional Riemannian manifold, (a, b run 
over 1,2,3. We restrict our discussion to three dimensional manifolds.) Our interest is in 
asymptotically flat spaces since we are interested in the energy and entropy of black holes. 
The definitions of energy, entropy and black holes all need an asymptotic region. The total 
energy or ADM mass of an initial data set is only well defined if an asymptotic structure 
(either flat or AdS) is fixed. Black holes are defined as regions of spacetime from which 
escape to infinity is impossible and thus refer to an asymptotic structure. We require that 
the metric tend to a fixed flat metric 6ab at infinity hab Sab + 0{l/r). Given an initial 
metric hab, the Ricci flow evolves the metric according to its Ricci tensor. The evolution 
parameter is r and the family of metrics on (S, hab{T)) satisfies the Ricci flow equation 

dha, 



Or 



-2Rab. (1) 



In the neighborhood of a point p E S, we can introduce a Riemann normal co-ordinate system 
and then the form of ([1]) becomes parabolic (V^ is the Laplacian in local co-ordinates) 

^ = ^'hab (2) 

and looks like a heat equation for the metric coefficients. However, in a general co-ordinate 
system, the PDE ([1]) is a degenerate parabolic equation, because of its diffeomorphism 
invariance. 

More generally, we will be interested in the Ricci flow modified by a diffeomorphism 

dh 

-2Rab + Da^b + DbU (3) 



dr 

where ^'^ is any vector field on S which vanishes at infinity. Calculationally, it is convenient 
to consider the two terms separately, defining a pure Ricci flow ([T]) and a pure diffeo 

dhab 



dr 



Da^b + DbU (4) 
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In the mathematical hterature on Ricci flows (which deals with compact spaces) the term 
Xhab is sometimes added to the RHS of ([1]) to define a "normalised Ricci flow". Such a term 
is inadmissible in the present physical context as it would rescale the metric at infinity and 
violate our asymptotic requirement that the metric tends to a fixed flat metric at infinity. 
In the standard initial value formulation 25(] of general relativity, the basic variables are 



the induced metric hab on a spatial slice S in the space-time manifold M, S G M and the 
extrinsic curvature kab of S. We also use the notation k = h°'^kab- These variables are 
subject to constraints: 

Dbik''^ - h^^k) = 87rf (5) 

and 

R+k'^- kabk""^ = levrp, (6) 

where j'' is the matter current and p is the matter density. The matter is required to satisfy 
"energy conditions", the dominant, weak or strong energy condition. 

A three dimensional manifold (S, hab) can be viewed as a time symmetric intitial data 
set for Einstein's equations. By "time symmetric" we mean that the extrinsic curvature kab 
of S has been set to zero, so that (S, hab) is totally geodesic. (This is similar to choosing 
initial data in classical mechanics so that all the momenta vanish. Dropping N particles from 
rest is an example.) With kab (and the matter current j") set to zero, the diffeomorphism 
constraint ([5]) is automatically satisfied and the Hamiltonian constraint reduces to R = IGirp. 
A physically important constraint on initial data for general relativity is that the data are 
subject to an energy condition. The dominant, weak and strong energy conditions all imply 
the local energy condition, which states that the local energy density is non-negative. This 
translates into the geometrical statement that the scalar curvature R of hab be non-negative. 
The Ricci flow has the appealing property that it preserves the non-negativity of scalar 
curvature 2J]. This follows from the "maximum principle" for the scalar curvature. Under 



Ricci flow(l3]), the scalar curvature satisfies the non linear heat type equation 

— = V^R + 2R'''Rab + C^R, (7) 
or 

which implies that ^ is positive at a minimum of R. Eq.Q provides us with a flow on 
the space of initial data to Einstein's equations which remains within physically allowed 
(non-negative scalar curvature) data. 
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III. GEOMETRIC QUANTITIES OF INTEREST 



Let 5" be a closed surface in S, jij, = 1,2) the pull back or induced metric on S, 
TZ the scalar curvature of {S,'~fij) and K the trace of its extrinsic curvature. We will be 
interested in the evolution of some geometric properties of S under the Ricci flow. These 
are the area of S, 



AS) 



dA 



and the Hawking mass of S 



(8) 



(9) 



Our interest in these quantities stems from their physical significance. The area of apparent 
horizons is related to the entropy of Black Holes and the Hawking Mass is related to the 
Energy. The Hawking Mass of a surface S can be physically interpreted as the mass contained 
within the surface 5*. While there are some problems with this interpretation (positivity is 
not always assured), the Hawking mass is an useful notion 2l|, |22j of quasilocal mass. It 
vanishes in the limit that 5* shrinks to a point and becomes the ADM energy for a round 
sphere at infinity. In fact, the supremum of M.{S) over S is the ADM mass as one can see 



from 



21 



22l |. Unlike the ADM energy, which is only well-defined for asymptotic spheres. 



the Hawking mass is defined for any closed surface S. Under the RF, a surface 5, which is 
initially asymptotic may shrink into the interior of S. For this reason Hawking mass is a 
more convenient object to study than the ADM mass. 

Rather than the Hawking mass, it is more convenient to deal with the related dimension- 
less quantity the "compactness" of S 



C{S) = [ dA{2TZ - 
Js 



(10) 



The quantity C{S) has been 



22l | in their approach to positive mass 



which is a combination of the Hawking mass and the area 

used to good effect by Geroch, Jang and Wald [! 
;heorem and the Penrose inequality. In fact their work forms the base for recent progress 
2o| on the Riemannian Penrose inequality. We will see that C{S) tends to zero as S tends 

to a round sphere of infinitesimal radius and also as S tends to an asymptotic round sphere. 
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IV. SPHERICAL SYMMETRY 



The Ricci flow is a tensor evolution equation and therefore commutes with diffeomor- 
phisms. It follows that isometry groups are preserved under the Ricci flow. One way to 
approach the Ricci flow is to start with symmetric situations so that the complexity of the 
flow is reduced. If we impose so much symmetry that the spaces of interest are homoge- 



neous, the RF becomes an ODE rather than a PDE. Such situations have been studied |24l|. 
However, this assumption is too restrictive from our present physical motivation. We would 
like to deal with asymptotically flat (or AdS ) spaces since energy is deflned with respect 
to an asymptotic structure. Homogeneous spaces which are asymptotically flat would be 
everywhere flat and not very interesting. Spherical symmetry has the advantage that there 
are non trivial asymptotically flat spaces. The Ricci flow reduces to a PDE with just two 
independent variables, the r and r co-ordinates. As we will see, spherically symmetric spaces 
will provide us with an analytical as well as numerical testing ground and pave the way for 
the general treatment. This special case is useful since one can develop a physical feel for the 
geometrical quantities of interest and easily produce physical examples and counterexamples 
as a guide to intuition. 

While setting up the spherically symmetric initial data set, we will work with two forms 
of the metric. Each of these has its use and its limitation. We will refer to them as "a-form" 
and "b-form". They correspond to different choices of co-ordinate gauge. 

• The a-form: In this form of the initial data set, the 3-metric is taken as 

ds'^ = a{r)dr^ + r'^i.dO'^ + sin^edcj)'^). (11) 

With this form of the metric we calculate the Ricci tensor and the scalar curvature: 
The nonzero components of the Ricci tensors are (a prime means differentiation with 
respect to r) 

E 

ra 

a'r 1 

Ree = + 1 , 

2a^ a 

R^^ = sin^9{Rge). 



The scalar curvature is 



2 2a' 



The a-form is useful for the study of the evolution of Hawking mass A4h- However, 
this form of the metric is not useful if there is an apparent horizon because in that 
case a(r) blows up at the apparent horizon. To study the evolution of the area of the 
apparent horizon under RF we use the "b-form" of the metric discussed next. 

The b-form: In this form of the initial data set, the 3- metric is taken as 



ds^ = dr' + h{r){de' + sin'ed(j)'). (12) 

With this metric we again calculate the Ricci tensor and the scalar curvature: The 
nonzero components of the Ricci tensor are 

6'2 - 2bb" 



The scalar curvature is 



Rrr ' 
b" 

Ree ~ ^ ^ ~2' 
Rct><t> = sin^9{Ree)- 



b'^- ^b{b" - 1) 
R 1^2 • (13) 



a-form and Hawking Mass: Let us start with the a-form f|TT]) of the metric and evaluate the 
Hawking mass functional ([9]) for 5* chosen to be a sphere r = constant. For this spherical 
topology, the first term in ( |T0|) gives IGvr by the Gauss-Bonnet theorem and 



Mh{S) := ^^[IQti- I K^d A]. (14) 



647r2 V Js 
Let fia be a unit normal to the surface S. The normalization 

h"'^hahb = a~^{r)nrhr = 1 (15) 

fixes 

hr = \fa (16) 

and so 

n'^ = fi,0,0). (17) 



/a 

The trace of the extrinsic curvature is 



K = Dafi'' = (18) 



r\ a 



so we have 



and 



K^dA 



A 



47ri/2 



167r 

a(r) 



1 - 



a(r) 



(19) 



(20) 



For flat space a(r) = 1 and we have Mh{S) = for round spheres in flat space as expected. 
If we take the example of the exterior Schwarzschild space 



a(r) 



2M 



(21) 



and then for any r > 2M, 



Mh{S) 



Trr^ 



40? 



1-1- 



2M 



M 



(22) 



For a general a- form metric we can write a(r) as I 1 — 



2M(r) 



. We find that [25 1 



MffiS) = M{r) 



and the "compactness" works out to 



C{r) = 327r 



M(r) 



(23) 



(24) 



In order to get a better physical feel for what these geometrical quantities mean, let 
us consider some typical distributions of matter. Let us choose the matter density p(r) = 
(l/167r)i?(r) positive and plot the functions p(r), C(r) and Muir). Figure([T]) displays the 
forms of these functions (in arbitrary units) for a spherical shell of matter. C(r) increases 
to a maximum value and then decreases to at infinity. For a matter distribution of two 
momentarily static shells of matter, the slightly more complex behaviour of C(r) is shown 
in Fig. ([2]). For a star C(r) attains its maximum near the surface of the star. From fl2^ . 
in the Newtonian limit C{r) is a constant times the dimensionless Newtonian potential, or 
the mass to radius ratio. Hence the name "compactness" is justified. 

Note that Mii{r) monotonically increases with r to attain its asymptotic value. This 
is due to the local energy condition, which implies positive scalar curvature i? > 0. This 
condition is conveniently stated in terms of the function M(r). The scalar curvature of 
( ITTj) is given by i? = 4M'(r)/r^ and so the constraint of positivity of scalar curvature 
simply states that M(r) is a non-decreasing function of r. Assuming that the form (fTTj) 



holds all the way to the origin and that the scalar curvature R is finite, we have M(0) = 
0,M'(0) = 0, M" = 0,M"'(0) > 0. M(r) increases from zero and tends to an asymptotic 
value Madm = linir->oo M(r) which is the ADM mass of the metric. In a sense 25|], M.H{f) 
measures the total mass contained within a sphere of areal radius r. M.H{r) is non-negative 
for all r and and from fl20p . we conclude that a(r) > 1 for all r. 

Note that from (fT9!) it follows that if the space contains an apparent horizon {K = 0), 
a(r) must diverge. This is why the a-form is unsuitable for treating apparent horizons. The 
b-form does not suffer from this problem. 

Geometric quantities in the b-form: We now consider a round sphere C S given by 
r = constant in the spherically symmetric "b-form" (IT^ of the metric. The unit normal 
satisfies 

h°'^hahb = firhr = 1 (25) 
and = (1, 0, 0). The area ^ of 5 is given as 

A{r)= [ ^ddd(P = A'Kh{r), (26) 
Js 

where 7 = h'^sin'^9 is the determinant of the induced metric on S. The trace of the 
extrinsic curvature is given by 

K := Dafi'' = p (27) 

(where a prime indicates differentiation with respect to r). The general formula for the 
compactness reduces in spherical symmetry to, 

C{r) = 167r - / y/^dOdcPK^ = 16n — (28) 

Js b 

and the Hawking mass is 

Note that Ain depends on b{r) and its derivative b'{r), in constrast to the simple algebraic 
relation fl20p we had in the a-form of the metric. 

Area under Ricci flow: We consider a metric initially in the b-form, evolving under a pure 
Ricci fiow (without a diffeomorphism term). Under this evolution, the b-form may not be 
preserved. We view S* as a fixed surface of S and so the co-ordinate location of the surface 
S does not change and hence ^ = 0. From the Ricci fiow we have 

dhf)o 



dr 

10 



-2Ro 



from which follows the formula 

dA dh ^ ,^ 

= 47r— = -47r(2 - b ) 

dr or 

for the instantaneous rate of change of area of 5*. 

Using the scalar curvature R for the "b-form" of the metric ( |T3l) we see that, in spherical 
symmetry, 

= f V^ded<pR-^C (30) 
and so we arrive at the inequality (since -R > 0) 

In the case of the Schwarzschild space, R = and so the first integral in (!30l) vanishes and 
we have 

Thus Schwarzschild space saturates our inequality (|3T1) . just as it saturates the Penrose 
inequahty. 

Area of apparent horizons under Ricci Flow: Let S* be a minimal surface (or apparent 
horizon, they coincide in the case of time symmetric data) in S z. e S* is a closed two manifold 
embedded in E with the property that the trace of the extrinsic curvature vanishes. We want 
to see how the area of S varies under the RF. We start with the spherically symmetric "b- 
form" of the metric (fT2|) . Let the location of the apparent horizon be at r = rg. From (127!) . 
we have that K = b' /b\rQ = The condition that the surface r = ro be an apparent horizon 
is 

b'\r=ro = 0. (33) 

Can a minimal surface spontaneously appear if none was present initially? The answer 
is no, as the following argument shows by contradiction. Let us consider the b form of the 
metric and evolve it by the Ricci flow supplemented by a suitable radial diffeomorphism i^" 
chosen to maintain the 6-form 

dh 

_nL = -2R„ + 2{Di)rr = 0, (34) 

^ = |^ = -2i?. + 2(DO... (35) 
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We may write = drf and using (IMI) and choose / so that 

^„ _ b'^ - 2bb" 



262 ■ 

(This leaves some freedom in /, but this does not affect the following.) (135!) then gives us 
the evolution equation for b{r, r) in the b form of the metric 

— = 6 -2 + 6/. 

Differentiating this equation we arrive at an evolution equation for b' 

^ = b'" + b"f + b'f", (36) 

which gives us a maximum principle for b'. Suppose that 6' > for all r < tq and that for the 
first time r = Tq, a minimal surface appears (6'(ro) = 0) at r = tq. We have f^(ro) < 0, since 
b' decreased to zero. On the other hand, since tq is a minimum of b', we have 6"(ro) = and 
b"'{ro) > 0. A glance at ( l36l) shows that ^ > 0, which is a contradiction. Thus a minimal 
surface cannot spontaneously appear under RF if it was not initially present. 

Regions where 6' < are called trapped regions. Since the RF is continuous, such regions 
evolve continuously. Trapped regions are bounded by minimal surfaces {b' = 0). From 
the last paragraph it is clear that trapped regions cannot spontaneously appear under RF 
since they are accompanied by minimal surfaces. However, trapped regions can continuously 
shrink to zero and disappear. When this happens, the two minimal surfaces that form the 
boundary of the trapped region merge and disappear. More descriptively, a minimal surface 
{b' = 0, b" > 0) merges with a maximal surface (6' = 0, b" < 0) and disappears. 

Setting aside such mergers, let us study how the area of minimal surfaces evolves under 
the pure RF (1). During the RF the metric changes and the location of the horizon may 
change and so ro = ro('r) where r is the parameter of the RF. Also the geometry of 5* will 
change. In principle both these effects could lead to change of area. The area is given by 
(1261) A{r) = 47r6(r). The total rate of change of area, therefore, is 

'db 



An— = A-K 
dr 



dr 



drQ ^ db 



(37) 

r=ro- 

The first term vanishes because of the apparent horizon condition fl33|) . The second term is 
evaluated by specialising eq. (j3Til to an apparent horizon. So the area A = Anb satisfies 

— < -47r. 

OT 
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This implies that the area of the horizon is decreasing at least linearly with r. Since the 
area was finite to begin with, we find that (if the horizon persists) b evaluated at the horizon 
goes to zero in a finite r. 

Next we see that as 6 ^ we approach a singularity. To show this we will suppose that 
R is finite and will arrive at a contradiction. As 6 — > the scalar curvature 

P 2(1 - b") 

is finite only if b" = 1. We Taylor expand in powers of r — rg about tq, the location of the 
apparent horizon 

b{r) = b{ro) + 6'(ro)(r - Tq) + h"{ro){r - rof + .... 

As 6(ro) 0, ^"('"o) 1 and 6'(ro) = due to the apparent horizon condition ( 133|) . we have 

1 



so the metric is 



b{r) = |(r-ro)^ 



ds'^ = dr^ + ^(r - rofide"^ + sin^ 



Shifting the r co-ordinate r ^ (r — Tq) gives the form 

ds^ = dr'^ + ^r^{d9^ + sin^ 
So the volume of a ball of radius r is 

V r^sin^O 1 / 47rr^ 



4 4V 3 

We then find from the expression for the volume of a ball of radius r (as r ^ 0) centered at 
point p that 

volume B{p,r) = ( -;r] \ ~ 777: -^(p)'"^ 



3 y V 30 

and so the scalar curvature blows up, R{p) ~ r~^, which is a contradiction to the assumption 
of finite R that we started with. So an apparent horizon which persists under RF results in 
a singularity in a finite amount of r parameter r < Tq, where Tq = A{0)/{A7t). 

Compactness under RF: We have previously derived the relation (l24l) for the compactness 
of a sphere r = const from which follows 

C(r) = 167r[l - 1/a]. (38) 
13 



A calculation reveals that under RF 

dC 2M'(r) 2M" , , 

dr r 

While the first term here is of definite sign, the second term is not. As a result, the rate 
of change of C with r is not monotonic. Note however that for Schwarzschild M(r) is a 
constant M and so dC/dr = 0. 

The Hawking mass is given by A^//(r) = rC{r)/ (327r) and to work out its rate of change 
^^^^ under the RF, we need to find dr/dr. This is easily read off from (l30l) and we find 

dr -M' C 



dr r 32'7rr 



(40) 



Putting these equations together we can work out "^"^^ and it turns out to be a linear 
combination of M' and M" . As a result, the Hawking mass is not monotonic along the pure 
Ricci flow in the a-form of the metric. It may be that adding a suitable diffeomorphism will 
result in monotonic behaviour for M.h- 

Maximum principle for compactness: However, we can make a diffeo invariant statement 
about the behaviour of compactness. We have already seen that if a minimal surface is 
initially present this may result in a finite r singularity. Let us suppose that no minimal 
surface is initially present. None develops under RF, as we saw earlier. We can therefore 
use the a-form of the metric. We have seen that C(r) starts from at r = 0, reaches a 
maximum (perhaps several local maxima) and then decays to zero as r — > oo. Let us focus 
on Cmax the value of C at its absolute maximum. Just as for apparent horizons, the value 
of Cmax is not affcctcd by a diffeomorphism: moving the surface does not affect the value 
of Cmax since we are at a maximum. We choose a diffeomorphism to preserve the a-form of 
the metric and then the evolution equation for a(r) reads 

da{r) _ a"{r) _ ?>{a'{r)f _ 2(a(r) - 1) + ra'(r)(l - a{r))/a{r) 
dr a{r) 2a{rY ' 

where a prime denotes differentiation with respect to r. We have studied the evolution of 
(14T|) using numerical techniques. Some of the results presented here were initially suggested 
by the numerical evidence. 

Let us focus on the maximum value of a(r). Recall that a(r) > L At the maximum value 
of a(r), we have a'{r)\max = and a"{r)\max < 0. So from equation fllTl) we see that the 
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maximum value of a(r) is monotone non-increasing as the flow parameter r increases i.e., 

'9a(r)^a^ _ a"{r) _ 2{a{r)max - 1) ^ ^ /^2) 
5r a(r) ~ 

A maximum of C corresponds via (138!) to a maximum of a(r). Clearly, a maximum principle 
for a(r) implies a maximum principle for compactness C{r) and we have the inequality 

dC{r) 



/max 



dr 



< 0. (43) 



We can use this maximum principle to comment on the long time existence of the spherically 
symmetric asymptotically flat Ricci flow. Let us suppose that no minimal surface is initially 
present. From fH2]) . we see that for the LHS of fH21) to vanish we must have a{r)max = 1, 
which implies a(r) = 1 identically. This describes flat space. If the initial metric is not 
flat space, its a{r)max decreases with the flow and flnally attains the flat space flxed point 
a(r) = 1. Our argument shows clearly that in spherical symmetry, the only asymptotically 
lat flxed point of the flow is flat space. This recovers results obtained by other methods 







To summarise, the maximum principle for compactness leads to a criterion for the exis- 
tence of the Ricci flow in the asymptotically flat case. If there are no apparent horizons, the 
flow exists for all r and converges to flat space. If there are apparent horizons, the Ricci 
flow either terminates in flnite time singularity or removes the horizons by mergers. 



V. AREA UNDER RICCI FLOW 



We now drop the assumption of spherical symmetry and deal with a general asymptot- 
ically flat manifold S, with one end at inflnity and a flxed closed orientable surface S of 
arbitrary topology embedded in S. The induced metric of 5* is written jij and extrinsic 
curvature tensor of S is written Kij, where i,j are two dimensional indices in the space 
tangential to S. We will sometimes use projected a, b indices for these. The trace of the 
extrinsic curvature is written K = Y^Kij. As the metric hab evolves according to the Ricci 
flow ([1]), how does the area change with r? Since the metric of S is changing, we have to 
remember that the unit normal to S is also changing with the metric. 

Let us deflne the surface S as the level set of a function 77 on E which is strictly increasing 
outward from S. Quite independent of any metric, the normal rja = DaTj is a well defined 
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CO- vector. r]a is non-zero (since we assumed t] is not locally constant). (Choosing a different 
function rj will result in multiplication of rja by a positive function on 5*.) The unit normal 

Va 



(44) 



(77.77)^/2 

depends on the metric. As the metric changes the unit normal can only change by a multiple 
of itself dfia/dt = aha- By differentiating h.h = 1, we arrive at 



Idk 



ab 



a 



2 dr 



-naUb, 



(45) 



(where in our notation, is defined as with its indices raised using the metric tensor 



h 



ab\ 



dA 



Starting from the definition (IHl) we compute 



— — = I d X 

dr 



d_V2 

dr 



We easily see that 



and we arrive at 



dr 2^^^ dr 



Area under ricci flow: Using the form ([T]) of the Ricci flow, we find 

^ = j^V^d'x[h-h'Rat-R]. 



From the Gauss-Codazzi equation 251] , we have 



n 



^h^'Rab -R= -1/2 [i? + n+ {KijK'^ - K% 



(46) 

(47) 
(48) 

(49) 
(50) 



which can be rearranged to give 



4 Js 



(51) 



The second integral in (15T!) is identified as —C{S)/A, minus one fourth the compactness 
integral of S and the first integral, which is of definite sign can be dropped to arrive at the 
inequality 

dr ~ A ^ ^ 
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This inequality is one of the main results of this paper. This result can be reexpressed in 
terms of the Hawking Mass: 

< ~24n'/^MH{S). (53) 

dr 

Thus the rate of decrease of area of a closed 2-surface under Ricci flow is bounded by the 
Hawking mass. As we mentioned earlier, the inequality ( l52l) is saturated in the case of 
the spheres of Schwarzschild space (which is given by ffTTl) with a(r) = (1 — 2M/r)~^). In 
this case R = and the spheres are shear free {Kij = ^Kjij), so the first integral in (15T|) 
vanishes. We conjecture that round spheres in Schwarzschild are the only surfaces which 
saturate this bound. 

As a simple application of this inequality, let us consider fiat space. Since the Ricci 
tensor vanishes we have that dA/dr = and so the LHS of ( |53l) vanishes. We arrive at 
the conclusion that for all surfaces in fiat space, the Hawking mass is non positive! This 



fact has also been noticed in 28|, where a direct proof is given. In fact, the converse of this 
statement is also true: Given positive scalar curvature, fiat space is the only one for which 
the Hawking mass is non-positive. To see this, note that the supremum over S" G S of the 
Hawking mass is the ADM mass and if this supremum vanishes, it follows from the positive 
mass theorem that the space must be fiat. 



Area under diffeos: Under a diffeo, the metric changes as in (jl]) and so we have from fH8|) 
^ = i y ^d^x{h-' - h-h^)2D,i, = j ^d^x[DaC - n'^n^'D^i,] (54) 
If we suppose that is tangent to then with Da denoting covariant derivative intrinsic 

to (7ab, S), 

Daib = i/i,''DAy (55) 



gives us 



So, as one would expect. 



DaC = r'Da^b = (h'^' - n-n')Daib (56) 



^d^xlDaC] = 0, (57) 



since this is a divergence over a closed surface S. It is therefore enough to consider the 
component of normal to S. Let 

C = un^. (58) 
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Then 

dA 

dr 



j ^d^xih"'' -n''n'')Da{uhb), (59) 



which works out to the standard answer 



j ^d^x u-i^'^DaUi, = j ^d^x uK. (60) 



dA 

Area of horizons under ricci flow: If S* is a minimal surface, we find from (fTOi) that since 
K = 0, the area changes according to 

dA{Sl < IciS) = (61) 
dr - 4 ^ ' 4 ' ^ ^ 

where x{S) is the Euler characteristic of S. For a minimal surface of spherical topology, we 
have 

^ < -4-. (62) 
which we had seen earlier for the special case of spherical symmetry. This result is unaffected 
by adding a diffeo to the RF because of (l60ll and the minimal surface condition K = 0. 

An interesting special case is one for which S is the outermost horizon with respect to 
asymptotic infinity. S is defined as the boundary of the region having trapped surfaces. S 
is a minimal surface and is known to have spherical topology. According to (!62l) . the area 
of S will shrink. Under the Ricci flow, the trapped region cannot disappear suddenly, but 
evolves continuously, because trapped surfaces remain trapped under small perturbations. 
Near S there will be a new minimal surface with the same area. If the initial area is it 
dissappears within a time Aq/Atc, either by merger or by shrinking to zero. 



VI. HAWKING MASS UNDER RICCI FLOW 

The Hawking mass of a closed surface S given by ([9]) is a combination of the "com- 
pactness" of 5* and its area. To study the evolution of the Hawking mass, it is enough to 
understand the evolution of the compactness, the evolution of area being already treated in 
the last section. From the formula ( JTOl) for compactness, we see that the first term j d'^x^/jTZ 
drops out on differentiation since it is a topological invariant by the Gauss-Bonnet theorem. 
The second term gives 

dC{S) df^.. f.r.^dK^, fr.2dV^, 



dr dr 



j K'^d'x = - j 2K^^d'x - j K'^d'x. (63) 
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Using (H7I) for the second term in (l63l) and the formula 

dK _ d 
dr 



.{Dan'') 



dr dr 
we find after a straightforward calculation that 



dr ' 



dC 
dr 



- / dA- 



dAK 

K2 V 



dh, 



ab 



dr 



dr 



dr 



ricUd 



(64) 
(65) 

dr - ■ " dr ■ 

Equations ( l66ll give the general evolution of the compactness under for any one parameter 
family hab{r) of metrics. We will specialise these equations to the two cases of interest: the 
pure Ricci flow and a pure diffeo. 

Compactness under ricci flow: Substituting = —2Rab into fl66l) we find after simpli- 
fication using the contracted Bianchi identity that 
dC f 

— = / dA{K\R + n'^n'Rab) - 2K[2R^'Dank - h''Da{R''''n,hd)]}. (67) 

s 

Compactness under diffeos: How does the compactness change under a pure diffeo gen- 
erated by As we saw for the area, the tangential component of does not cause any 
change in the integral. The diffeo vector field can be characterised by its normal component 
£f_= uh°'. Substituting (j4]) into ( !66l) and simplifying and using the Gauss-Codazzi equation 



25| 



- 2nVRab + R = n+ {k'^ki 



(68) 



gives us the formula 22 1 



dC f 1 

— = / [2Kb''baU + uKa'^aij + ukR - -uK{2n - K^)]J^d^x 
Is 



(69) 



dr Js' ' """"" ' 2 

where Da denotes the intrinsic covariant derivative operator within the surface and cr*-' = 
_ If there exists a diffeo such that uK = 1 (this is the inverse mean curvature 

flow (IMC)) one can conclude [22] that under this diffeo 

(70) 



This pequality (which too ,s saturated by the Schwa.sehi,d space) was used by Jang and 

Wald21^ to prove the Penrose inequality. We can get an expression for the evolution of 
Hawking mass. The Hawking mass is 



Mh{S) 



647r3/2 



■C{S). 



(71) 
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We then have 

Knowing ^ and ^ for a flow, can calculate -^Mh- (l60i) and Geroch's inequality ( 1701) then 
imply that the Hawking mass is monotonic under the IMC flow. 

A maximum principle for compactness? The compactness C{S) of a closed surface 5* 
tends to zero for S tending to a small round sphere, and also tends to zero for a round 
spheres at asymptotic infinity. We would expect based on our experience with spherical 
symmetry, that somewhere in between, there is a surface 5* for which C{S) attains its global 
maximum. If this maximum is equal to IGvr, then 5 is a minimal surface and conversely. If 
the maximal compactness is less than IGtt, we may expect from the spherically symmetric 
case that this value C{S)\max "will monotonically decrease to zero. 

< 0, (73) 

ctr 

with the equality holding only for flat space. A zero value for C{S)\rnax would also imply a 
non-positive value for M.h{S). Supremising this over S tells us that the ADM mass of the 
space vanishes, which implies (by the positive energy theorem 29l|) that space is flat. We 
conjecture that fl73l) is true generally, but have not been able to establish the truth of this 
conjecture. If this conjecture is true, it would imply that in the long time, an initial metric 
without minimal surfaces would approach flat space. 



VII. CONCLUSION 



We have described some applications of Ricci flow techniques to asymptotically flat spaces 
in general relativity. Our main result is that under the Ricci flow, the rate of change of area 
of a closed surface is bounded by its Hawking mass. This inequality is saturated by the 
Schwarzshild space. Since the Schwarzschild space saturates the Penrose inequality as well 
as Geroch's fl70|) . our inequality may be related to these. We have also studied the behaviour 
of compactness under the Ricci flow. Our work in spherical symmetry suggests that there 
may be a maximum principle for the compactness. The compactness is a functional on 
closed two dimensional surfaces. The conjecture is that for positive scalar curvature, the 
maximum value of this functional decreases under the Ricci flow plus diffeomorphisms. This 
proposal is a functional maximum principle unlike the more usual ones which are which are 
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formulated on functions. If such a principle does exist, it would be of great interest as a 
diagnostic for the long term existence of the asymptotically fiat Ricci flow. Most of the 
mathematical work in the area of Ricci flows is concerned with closed three manifolds. One 
of our main points here is that there may be interesting physical applications of these ideas 
to asymptotically fiat spaces. For instance, |^ poses the question of stability of Euclidean 
fiat space under RF. If a maximum principle for compactness exists, it would imply that 
flat space is indeed stable under RF. 

In general a minimal surface can have any topology. Special interest attaches to the 
case where S is an outermost horizon. In this case one can physically identify S as the 
boundary of a black hole region. Outermost horizons have the property that the surface has 
minimum area (and not just stationary area). This immediately implies js^ (when curvature 
is positive) that the outermost horizon has spherical topology. Our general result shows 
that outermost horizons always shrink under the Ricci flow. This result is diffeomorphism 
invariant. If one identifies the area of an outermost horizon as black hole entropy, we 
arrive at the conclusion that entropy is monotonically decreasing along the RF. A similar 
conclusion has been reached by in a slightly different context: entanglement entropy of 
a two dimensional black hole is monotonically decreasing along the RF. 

If one uses the Ricci flow to model the approach of a system to thermal equilibrium, 
one would expect that the entropy increases along the flow. However, we have seen that 
outermost horizons have spherical topology and their area decreases under the flow. In fact, 
it has been observed 12|] that while the Ricci flow does have something in common (memory 
loss) with approach to equilibrium, the analogy is not perfect. A slight modification 12[ of 
the Ricci flow is necessary in order for the entropy to be identified with black hole entropy. It 
may be that such a modification of the RF is necessary for application to black hole physics. 

This paper was concerned exclusively with the case of asymptotically flat spaces. It 
would be interesting to generalise the treatment to allow for asymptotically AdS spaces. 
This appears to be a straightforward generalisation and some of the necessary changes are 



mentioned 



in Q, 



3Q|. 
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FIG. 1: Note that Mnif) increases with r to attain its asymptotic ADM value. But C(r) increases 
to a maximum value and then decreases to at infinity. 
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FIG. 2: Two shells of matter. 
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